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We present measurements of the CP-violation parameters S and C for the radiative decay 

-qKsj; for B rjKj wo also measure the branching fractions and for B^ r/K'^'y the time- 
integrated charge asymmetry Ach- The data, collected with the BABAR detector at the Stanford 
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Linear Accelerator Center, represent 465 x 10^ BB pairs produced in e'^e annihilation. The re- 
sults are S = -O.lSto H ± 0.12, C = -0.32t°;*° ± 0.07, 6(3° riK^j) = {l.lttl ± 0.4) x 10-^ 
B{B+ -> ■qK+'y) = (7.7 ± 1.0 ± 0.4) x 10~^ and Ah = (-9.0t9''g'' ± 1.4) x 10"^ The first error 
quoted is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Radiative B meson decays have long been recognized 
as a sensitive probe to test the standard model (SM) and 
to look for new physics (NP) [11,0. In the SM, flavor- 
changing neutral current processes, such as 5 — > 57, pro- 
ceed via radiative loop diagrams. The loop diagrams 
may also contain new heavy particles, and therefore are 
sensitive to NP. The measured branching fractions of in- 
clusive 6 — > 57 and exclusive radiative B decays are in 
agreement with SM predictions 0, 0, Recent esti- 
mates of the branching fraction of the inclusive 5 — > S7 
decay are affected by a theoretical uncertainties as large 
as the experimental ones 

In the SM the photon polarization in radiative decays 
is dominantly left (right) handed for b (6) decays, result- 
ing in the suppression of mixing-induced 0] CP asym- 
metries. Because we do not measure the photon helicity, 
we sum the decay rates for the left-handed and right- 
handed helicity states. Observation of significant CP- 
violation in these radiative decay modes would provide 
a clear sign of NP We search also for direct CP 
asymmetry in charged B decays, measuring the charge 
asymmetry Ach = (P" - P+)/(r" -I- P+), where P is the 
partial decay width of the B meson, and the superscript 
corresponds to its charge. This asymmetry in the SM is 
expected to be very small Q . 

In this letter, we present the first measurement of 
the mixing-induced CP violation in the decay mode 
B^ T]K'^j. Branching fractions for the decay modes 
B^ TjK'^j and 5+ rjK^j Q and time-integrated 
charge-asymmetry for B^ — > r]K~^j have been measured 
previously by the Belle and BaBAR 11 | Collabora- 
tions. We update our previous measurements with a data 
sample that is twice as large. 

The results presented here are based on data collected 
with the BaBAR detector 12| at the PEP-II asymmetric- 
energy e^e~ collider [l3| located at the Stanford Linear 
Accelerator Center. We use an integrated luminosity of 
423 fb~^, corresponding to 465 ± 5 million BB pairs, 
recorded at the T{AS) resonance (at a center-of-mass en- 
ergy of Vs = 10.58 GeV). 

Description of the BaBAR detector and of the recon- 
struction of charged and neutral particles can be found 
elsewhere [14| . The B decay daughter candidates are re- 
constructed through their decays tt" 77, rj ^ (^77)7 



and rj 



(Vstt)- Reconstruction and selection 



criteria of charged and neutral mesons, and primary pho- 
tons and the study of continuum and BB backgrounds 
are described in our previous paper [llj] . 

A B meson candidate is reconstructed by com- 



bining an r] candidate, a charged or neutral kaon 
and a primary photon candidate. It is characterized 
kincmatically by the energy-substituted mass toes = 
^ (s/2 -I- pg • pb)'^/Eq — and energy difference AE = 
Eg — ^^/s, where the subscripts and B refer to the ini- 
tial T(4S') and to the B candidate in the lab-frame, re- 
spectively, and the asterisk denotes the T{AS) rest frame. 
We require 5.25 < toes < 5.29 GeV/c^ and |A£;| < 0.2 
GeV. _ 

From a candidate BB pair we reconstruct a B'^ de- 
caying into rjK^"i {Brec)- We also reconstruct the decay 
point of the other B meson (i?tag) and identify its fla- 
vor. The difference At = tree — ^tag of the proper de- 
cay times tioc and itag of the reconstructed and tag B 
mesons, respectively, is obtained from the measured dis- 
tance between the Brec and Stag decay vertices and from 



the boost 
tribution [L 



F{At) 



?7 = 0.56) of the e+e 
I is given by: 



system. The At dis- 



\At\/r 



At 



-[1 =F Aw± 



(1 - 2w) {Ssm{AmdAt) - C cos{AmdAt))]. (1) 

The upper (lower) sign denotes a decay accompanied by 
a B^{B'^) tag, T is the mean B^ lifetime, Ato^ is the 
mixing frequency, and the mistag parameters w and Aw 
are the average and difference, respectively, of the prob- 
abilities that a true B^ is incorrectly tagged as a B^ 
or vice versa. In the flavor tagging algorithm [l^ there 
are six mutually exclusive tagging categories of different 
response purities and untagged events with no tagging 
informations. Tagging and At informations are used for 
the measurement of the CP- violation parameters S and 
C in the decay mode i?" rjK^'y. 

We use the same technique developed for B^ n^K^j 
decays [15[ to reconstruct the B° rj^^K'^^ decay point, 
using the knowledge of the K'^ trajectory and the average 
interaction point in a geometric fit. The extraction of At 
has been extensively validated in data [l3| and in a full 
detector simulation. In about 70% of the selected events 
the At resolution is sufficient for the time-dependent CP- 
violation measurement. For the remaining events the At 
information is not used. For both rj^-yK'^j and 773^ ^^^7 
modes we use the events which satisfy the requirements 
I At I < 20 ps and a At < 2.5 ps, where a At is the per-event 
error on At. 

We obtain signal event yields and CP-violation param- 
eters from unbinned extended maximum-likelihood (ML) 
fits. We indicate with j the species of event: signal, qq 
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continuum background, BB peaking background {BP), 
and BB non-peaking background (BNP). The input ob- 
servables are uies, ^E, the output of a Neural Network 
{NN), the r] invariant mass m^, and At. The NN com- 
bines four variables: the absolute values of the cosines of 
the polar angles with respect to the beam axis in the 
T(4S') frame of the B candidate momentum and the 
B thrust axis, the ratio of the second and zeroth Fox- 
Wolfram moments and the absolute value of the 
cosine of the angle 9t between the thrust axis of the B 
candidate and that of the rest of the tracks and neutral 
clusters in the event, calculated in the T(4S') frame. 

For each species j and tagging category c, we define a 
total probability density function (PDF) for event i as 



3,c 



r,imBsi-V,iAE^)-r,iNN^) 



(2) 



The factored form of the PDF is a good approximation 
since correlations between input observables are small. 
With Tij defined to be the number of events of the species 
j and fj^c the fraction of events of species j for each 
category c, we write the extended likelihood function for 
all events belonging to category c as 



Cc = exp ( - 



No 



JJ^ ('T'sig /sig, c'Pgig , c 



'^qqfqq,c'Plq + nBNP JbNP^cPbNP 

-nspfBP.c'PBp): 



(3) 



where rij^c is the yield of events of species j found by 
the fitter in category c and Nc the number of events 
of category c in the sample. We fix /sig,c, /bnp^c, and 
/bp,c to fBciav,c, the values measured with a large sam- 
ple of i3-decays to fully reconstructed fiavor eigenstates 
(Sflav) [l9^|. The total likelihood function Cd for decay 
mode d is given as the product over the seven tagging cat- 
egories. Finally, when combining decay modes we form 
the grand likelihood C = Il'^c'- 

The PDF Vsig{At, (JAt',c), for each category c, is the 
convolution of F{At; c) (Eq. [T|) with the signal res- 
olution function (sum of three Gaussians) determined 
from the -Bfiav sample. The other PDF forms are: the 
sum of two Gaussians for 7'sig(™Es)j ^sig(A£^), and 
Psigirfiri)', the sum of three Gaussians for Vqq (At), 
'PsNpiAt), and 7^Bp(At); a non-parametric step func- 
tion for VjiNN) [iOl; a linear dependence for Vqq {AE), 
Vbnp{AE), and Vbp{AE); a first-order polynomial plus 
a Gaussian for Vqq {rriji), VBNpifnri)-, and VBp(TTiri)] and 
for Vqq (ttt-es), VBNp{mEs), and VspimEs), the func- 
tion a;-\/l — exp [—^(1 — x"^)] , with x = 2mEs/^/s PH . 
where for the BP PDFs we add a Gaussian. We allow 
qq background PDF parameters to vary in the fit. 

We determine the PDF parameters from Monte Carlo 
(MC) simulation for the signal and BB backgrounds. 



while using sideband data (5.25 < mss < 5.27 GeV/c^; 
0.1 < |A£:| < 0.2 GeV) to model the PDFs of contin- 
uum background. Large control samples of B decays to 
charmed final states with similar topology and a smear- 
ing procedure applied to photons during the event recon- 
struction are used to verify the simulated resolutions in 
niES and AE. Where the control data samples reveal dif- 
ferences from the MC in mass resolution, we shift or scale 
the resolution used in the likelihood fits. The largest shift 
in TOES is 0.6 MeV/c^. Any bias in the fit is determined 
from a large set of simulated experiments in which the 
qq and BNP backgrounds are generated from the PDFs, 
and into which we have embedded the expected num- 
ber of BP and signal events chosen randomly from fully 
simulated MC samples. 

We compute the branching fractions and charge asym- 
metry from fits made without At or fiavor tagging. The 
free parameters in the fit arc: the signal, qq , BNP 
and BP background yields; the bin weights of the step 
function for Vnq {NN); the slopes of Vqn (AE) and 



qq 



); ^; and for charged modes the signal and back- 
ground Ach- We apply the same procedure to extract S 
and C. In this case we add in the fit the At variable 
and the flavor-tagging information. As free parameters 
we have also S, C, and the parameters of the Vqq {At) 
PDF. 

TableUlists the results of the fits. The corrected signal 
yield is the fitted yield minus the fit bias which is in the 
range 2 — 4%. The efficiency is calculated as the ratio 
of the number of signal MC events entering the ML fit 
to the total generated. We compute the branching frac- 
tions from the corrected signal yields, reconstruction ef- 
ficiencies, daughter branching fractions, and the number 
of produced B mesons. We assume that the branching 
fractions of the T(4S') to B^B~ and B^B^ are each equal 
to 50%. We combine results from different channels by 
adding their likelihood functions, taking into account the 
correlated and uncorrelated systematic errors. 

The statistical error on the signal yield, S, C and the 
signal charge asymmetry is taken as the change in the 
central value when the quantity — 2 ln£ increases by one 
unit from its minimum value. The significance S {a) is 
the square root of the difference between the value of 
— 21n£ (with systematic uncertainties included) for zero 
signal and the value at its minimum. 

Figure [T] shows, as representative fits, the projections 
onto Toes and AE while Fig.[2]shows the projections onto 
At and the raw asymmetry between B'^ and tags. In 
these projections a subset of the data is used for which the 
signal likelihood (computed without the variable plotted) 
exceeds a threshold that optimizes the sensitivity. 

Figure [3] shows the distribution of the r]K invariant 
mass for signal events obtained by the event-weighting 
technique (sPlot) described in Ref. [22|. We use the co- 
variance matrix and PDFs from the ML fit to determine 
a probability for each signal event. The resulting distri- 
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TABLE I: Number of events in the sample, corrected signal yield, detection efficiency e, daughter branching fraction 
product Yl Bi , significance S (a) (including systematic uncertainties) , and measured branching fraction B with statistical error 
for each decay mode. For the combined measurements we give S (a) and the branching fraction with statistical and systematic 
uncertainty. For the neutral mode we give the S and C parameters for each decay mode and for their combination. For the 
charged modes we also give the measured signal charge asymmetry Ach- 



Mode 


N 


Yield 


e (%) UB^{%) S (a) 


s(io-^) 


Ach (10-2) 


S 


C 




3690 
2282 




12 
10 


13.6 
7.8 


3.3 
2.1 
3.9 


r.itll ± 0.4 




-0.04 ±0.62 
-0.45 ±0.81 


-0.24 ± 0.44 
-0.71 ±0.87 
-0.32^039 ± 0.07 




11620 
10738 




19 
14 


39.4 
22.4 


6.5 
4.5 
8.0 


'•°-i.o 
7 4+1.7 

'•^-1.6 

7.7 ± 1.0 ± 0.4 


-4± 12 
-24 ±20 
-9.01^",^ ± 1.4 
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FIG. 1: The B candidate rriES and AE projections (see text) 
for rjK'^'y (a, b), r/K^j (c, d). Points with error bars (statisti- 
cal only) represent the data, the solid line the full fit function, 
and the dashed line its background component. 




At (ps) 



FIG. 2: Projections (see text) onto At of the data (points with 
error bars), fit function (solid line), and background function 
(dashed line), for (a) and (b) tagged events, and (c) 
the raw asymmetry [Ngo — N-go ) / {Ngo + N-go ) between i?" 
and B° tags. 



butions (points with errors) are normalized to the signal 
yield. This mass distribution is useful to compare with 
theoretical predictions for radiative decays. 

The main sources of systematic uncertainties for the 
time-dependent measurements come from the variation 
of the signal PDF shape parameters within their errors 
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FIG. 3: Plot of rjK invariant mass for signal using the weight- 
ing technique described in the text for the combined sub-decay 
modes: (a) B+ vK^l, (b) B° r/K""/. Errors are statis- 
tical only. 



(0.08 for S, 0.04 for C), and from BB backgrounds (0.09 
for 5, 0.06 for C). Other minor sources are SVT align- 
ment, beam spot position and size, and interference be- 
tween the CKM-suppressed b —^ ucd amplitude and the 
favored b —^ cud amplitude for some tag-side B decays 
[m . The Bfiav sample is used to determine the errors as- 
sociated with the signal At resolutions, tagging efficien- 
cies, and mistag rates. Published measurements [4j for 
T and Arrid are used to determine the errors associated 
with them. Summing all systematic errors in quadrature, 
we obtain ±0.12 for S and ±0.07 for C. 

The main sources of systematic uncertainties for the 
branching fraction measurements include uncertainties in 
the PDF parameterization and ML fit bias. For the sig- 
nal, the uncertainties in PDF parameters are estimated 
by comparing MC and data in control samples. Varying 
the signal PDF parameters within these errors, we es- 
timate yield uncertainties of 3-23 events, depending on 
the mode. The uncertainty (1-3 events) from fit bias 
is taken as half the correction itself. Systematic uncer- 
tainties due to lack of knowledge of the primary photon 
spectrum are estimated to be in the range 2-3% depend- 
ing on the decay mode. Uncertainties in our knowledge 
of the efficiency, found from auxiliary studies 
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elude 0.4% X Nt and 1.8% x N^, where Nt and N. 



m- 
are 



the numbers of tracks and photons, respectively, in the B 
candidate. There is a systematic error of 2.1% in the effi- 
ciency of K'i reconstruction. The uncertainty in the total 
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number of BB pairs in the data sample is 1.1%. Pub- 
lished data Q provide the uncertainties in the B daugh- 
ter branching fraction products (0.7-1.8%). 

A systematic uncertainty of 0.014 is assigned to Ach- 
This uncertainty is estimated from studies with signal 
MC events and data control samples and from calcula- 
tion of the asymmetry due to particles interacting in the 
detector. 

In conclusion, we measure the time-dependent CP vi- 
olation parameters in the decay mode B'^ — )■ rjK'^^: 
S = -0.18toil±0.l2 and C = -0.32t[Jj° ± 0.07. These 
results are consistent with no CP-violation in this mode. 
We also measure the branching fractions, in units of 10~^, 
B{B° r)K°-f)= 7.ltl[l ± 0.4 and B{B+ r]K+j) = 
7.7 ± 1.0 ± 0.4, in agreement with the results from 
Belle [lo| and the previous BABAR results [llj. The 
measured charge asymmetry in the decay B~^ —^ rjK^^ 
is consistent with zero. Its confidence interval at 90% 
confidence level is [—0.25,0.08]. 
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